System for Monitoring and Indicating Filter Life

ABSTRACT

Disclosed are systems, methods, and algorithms for monitoring and indicating filter life. In particular, the disclosed systems, methods, and algorithms may be utilized for monitoring and indicating the useful life of a filter in an internal combustion engine.

CROSS-REFERENCE TO RELATED APPLICATIONS

The present application claims the benefit of priority under 35 U.S.C. §119(e) to U.S. Provisional Patent Application No. 61/352,852, filed on Jun. 9, 2010, the content of which is incorporated by reference in its entirety.

BACKGROUND

The field of the invention relates to filters. In particular, the field of the invention relates to systems for monitoring and indicating filter life for a filter in an internal combustion engine.

SUMMARY

Disclosed are systems for monitoring and indicating filter life. In particular, the disclosed systems may be utilized for monitoring and indicating the useful life of a filter in a machine such as an internal combustion engine.

The disclosed systems utilize methods for monitoring and indicating filter life. The disclosed methods may be utilized for determining and indicating useful life of a filter for filtering fluid in a given machine. Typically, the methods optimize cost-efficiency and value to a customer or purchaser of a filter by maximizing filter service intervals and full-life usage of a filter and minimizing available but unused filtration capacity of the filter prior to the end of filter-life. The methods may maintain fuel economy by reducing or limiting the total work required to pump fluid through the filter. The methods utilize an algorithm for calculating remaining useful filter life as disclosed herein. Initial constants for the algorithm are based on laboratory performance of the filter, historical performance of the filter, or both. Further, in the methods, the initial constants for the algorithm may be adjusted during filter usage based on actual performance of the filter, behavior or use of the machine in which the filter is utilized, or both. The algorithm may be utilized to calculate remaining useful filter life based on parameters that include, but are not limited to ΔP, filter efficiency, and particle concentration downstream of the filter.

Suitable machines for the disclosed methods may include, but are not limited to, internal combustion engines and hydraulic machines. Suitable filters may include filters utilized for fuel-water separation, fuel filtration, lube filtration, coolants, hydraulics, and air.

Suitable machines for the disclosed methods may include an electronic control module (ECM) or a separate control module. In the methods, the initial constants for the utilized algorithm may be adjusted during operation based on behavior of the filter and resulting sensor input to the ECM which may be normalized or otherwise adjusted by the ECM, including historical behavior of the filter. In the methods, data may be obtained from the ECM and the obtained data may be used to calculate the volume of fluid flow through the filter. The volume of fluid flow thus calculated may be inputted into the algorithm utilized in the methods to indicate whether filter service is required. In some embodiments, the ECM may contain the algorithms and/or logic to do the analyses. In the disclosed methods, the initial constants for the algorithm typically may be adjusted during operation based on actual performance of the filter, behavior of the machine in which the filter is utilized, or both. The constants may be adjusted manually, automatically, or both. For example, in application, during operation, or in service, the initial constants for the algorithm typically may be adjusted during operation based on actual performance of the filter, behavior of the machine in which the filter is utilized, or both.

The initial constants may be adaptively varied or adjusted during operation, for example, if the observed pressure drop differs from the expected pressure drop based on the initial constants. For example, the initial constants may be adaptively varied or adjusted during operation based on parameters including, but not limited to, contaminant concentration, contaminant type, filter type, or a combination thereof.

In the disclosed methods, the initial constants may be adjusted during operation based on observed filter performance. Further, the algorithm utilized in the methods may be adaptively updated based on observed filter performance.

The disclosed methods may be utilized for determining and indicating useful life of a filter for filtering fluid in a given machine. In some embodiments, the methods provide an indication alert to a user that the filter is approaching an end-of-useful-life condition. For example, the methods may provide an indication alert to a user that the filter is approaching an end-of-useful-life condition based on total volume of fluid filter. Alternatively or in addition, the methods may provide an indication alert to a user that the filter is approaching an end-of-useful-life condition based on operating conditions to which the filter has been exposed.

In the methods the initial constants for the algorithm may be adjusted during operation during filter usage based on behavior or use of the machine in which the filter is utilized. For example, the initial constants may be adjusted during operation according to the number of start-up and shut-down events of the machine.

The algorithm in the methods may be adaptively adjusted for scheduling filter changes. In some embodiments, the algorithm is adaptively adjusted such that the filter is not replaced too soon, which may increase cost of operation for a machine that utilizes the filter. In other embodiments, the algorithm is adaptively adjusted such that the filter is not replaced too late, which may result in possible damage to the machine.

In the methods, the algorithm may be utilized as an adaptive tool which may be adjusted to more accurately estimate remaining useful life of said filter. The algorithm may be adjusted manually by service personnel. Alternatively or in addition, the algorithm may be adjusted automatically based on field experience to more accurately estimate remaining useful life of said filter.

The methods may include determining the volume of fluid filtered by the filter and using the algorithm to calculate remaining useful life of the filter based on the volume of fluid filtered. For example, the volume of fluid filtered may be calculated by tracking flow rate of fluid through the filter and by tracking time.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1 illustrates a schematic diagram of contamination sources that provide bases for one embodiment of an algorithm utilized in the systems and methods for monitoring and indicating filter life as contemplated herein.

FIG. 2 provides a graphic illustration of the effect of variable “A” on the behavior of ΔP.

FIG. 3 provides graphic illustrations of the effect of variable “B” on the behavior of ΔP. A. Linear plot. B. Logarithmic plot.

FIG. 4 provides graphic illustrations of the effect of variable “C” on the behavior of ΔP. A. Linear plot. B. Logarithmic plot.

FIG. 5 illustrates plugging relationships based on equations ΔP=A+B*exp(C*M) and ΔP=A+B*M^(C).

FIG. 6 illustrates a schematic diagram of one embodiment of an algorithm utilized in a system for monitoring and indicating filter life as contemplated herein including correction steps.

FIG. 7 illustrates a schematic diagram of one embodiment of an algorithm utilized in a system for monitoring and indicating filter life as contemplated herein including automatic correction steps.

FIG. 8 illustrates a schematic diagram of one embodiment of an algorithm utilized in a system for monitoring and indicating filter life as contemplated herein including automatic correction steps which calculates the remaining fluid volume that can be filtered before its contaminant loading capacity is reached.

FIG. 9 illustrates embodiments of a system for monitoring and indicating filter life as contemplated herein.

DETAILED DESCRIPTION

The present systems, methods, and algorithms may be described herein using several definitions, as set forth below and throughout the application.

Unless otherwise specified or indicated by context, the terms “a”, “an”, and “the” mean “one or more.” For example, “a system” or “an algorithm” should be interpreted to mean “one or more systems” or “one or more algorithms.”

As used herein, “about,” “approximately,” “substantially,” and “significantly” will be understood by persons of ordinary skill in the art and will vary to some extent on the context in which they are used. If there are uses of these terms which are not clear to persons of ordinary skill in the art given the context in which they are used, “about” and “approximately” will mean plus or minus ≦10% of the particular term and “substantially” and “significantly” will mean plus or minus >10% of the particular term.

As used herein, the terms “include” and “including” have the same meaning as the terms “comprise” and “comprising.” For example, a “system including a sensor” should be interpreted to mean “a system comprising a sensor.”

As used herein, the term “the” may be used interchangeably with the term “said.” For example, “the system” may be interchangeably referred to as “said system.”

As used herein, a “user” may include someone who is using a filter (i.e., a “customer” or an “operator”). A “user” further may include someone is monitoring the useful life of a filter to determine when the filter should be replaced (e.g., “service personnel”).

The systems and methods disclosed herein may be utilize to monitor and indicate when the useful life of a filter is approaching its end based on the total volume of fluid that the filter has processed and/or the operating conditions to which the filter has been exposed. When the useful life of a filter is approaching its end, a user is alerted so that the filter may be replaced before the useful life of the filter has ended. The present systems and methods may provide advantages over systems and methods that utilize a change in pressure to determine when the useful life of a filter is approaching (i.e., “a ΔP-useful life indication”). In particular, the disclosed systems and methods may be utilized in fuel-water separators, fuel filtration, lube filtration, coolants, hydraulics, and air, in which filtration media effectiveness is a function of exposure to the total volume of fuel flowed through the filtration media regardless of particulate contaminant levels in the fuel, which is one parameter that can trigger a ΔP-useful life indication. In addition, under certain operating conditions, for instance significant number of start-up/shut-down events, the ability of filtration media to retain removed particles may be affected prior to a ΔP-useful life indication. In some embodiments, “a ΔP-useful life indication” may be reached when ΔP is greater than about 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, or 30 kPa.

The disclosed systems and methods may be utilized to achieve a number of goals. In some embodiments, the disclosed systems and methods may be utilized to inform a user of the condition of a filter at an earlier point than indicated by other indicators, such as low power for a fuel filter or a ΔP indicator light. Also, by utilizing the disclosed systems and methods, a user may be alerted to when a filter should be replaced based on an approaching end to its useful life rather than based on time or mileage of usage for the filter. Further, the disclosed systems and methods permit users to maximize a filter product's useful life and minimize cost. In further embodiments, the disclosed systems and methods may be utilized to prevent and/or discourage the use of a filter product beyond the product's useful life, which may potentially compromise equipment utilizing the filter product in a harsh environment or in other applications. The disclosed systems and methods provide a user (e.g., an operator of a machine that incorporates the filter or service personnel that replace the filter) with an adaptive tool for scheduling filter changes, such that filter changes are not scheduled too soon thereby increasing cost of operation, and such that filter changes are not scheduled too late thereby resulting in damage to a machine that incorporates the filter (e.g., an engine). In further embodiments, the algorithm utilized in the disclosed systems and methods may be adjusted by users (e.g., service personnel that replace the filter) or automatically based on field experience to more accurately estimate remaining useful life of filter.

The disclosed systems and methods may include or utilize an electronic control module (ECM) and data captured by an ECM to calculate the volume of flow through a filter based on algorithms disclosed herein. For example, the disclosed systems and methods may use data captured by an ECM to calculate the volume of flow through a filter based on algorithm disclosed herein in order to indicate whether the filter is ending its useful life and should be replaced. The disclosed algorithms may be adjusted based on chosen factors (e.g., engine start-up/shut-down events) in order to better indicate when a filter is ending its useful life and should be replaced. The disclosed algorithms also may be utilized to indicate to a user that certain engine performance parameters should be adjusted, for example, engine performance parameters that typically alert a user that a filter should be replaced such as power limit.

The presently disclosed algorithms may be used to calculate remaining useful life of filter based on fluid volume filtered. Flow rate through the filter and time may be tracked and used to calculate volume filtered. The algorithm constants may be adjusted (e.g., manually or automatically) based on field experience, in order to improve accuracy over time with respect to estimating when a filter should be replaced. In this manner, the disclosed algorithm provides a tool that allows a user (e.g., an operator or service personnel) to schedule filter service at more convenient times and, with use, improves accuracy of prediction with respect to end of useful life for the filter, thereby reducing costs and improving reliability. In manual operation, a user returns to a service center for appropriate algorithm adjustments, while in automatic operation, supporting sensors typically are present in the systems (e.g., on an engine) and are connected to a controller.

The disclosed algorithms may operate as follows. Initially, filter-life is defined in terms of miles or hours based on default values entered into the algorithm. These default values may be obtained based on lab testing or prior experience in order to estimate remaining filter or component life. In some embodiments, the algorithm constants may be adjusted automatically as the filter is used. For example, if the observed ΔP or duty cycle is not as expected, then algorithm constants can be modified accordingly. In other embodiments, the algorithm constants may be adjusted manually based on observed conditions. For example, if engine behavior is not as expected, or if upon change-out, a filter appears to have remaining useful life or appears to have exceeded useful life.

As the filter is used, flow through the filter and time may be measured, the volume filtered may be calculated, and the remaining useful life of the filter may be determined using the algorithms disclosed herein. The constant values utilized in the algorithms may be adjusted such that the system continuously refines constant values (e.g., manually or automatically, as previously discussed), such that the algorithm becomes more accurate over time and becomes tailored to specific machines or applications.

The algorithm assumes that removal efficiency for any given filter is a function of ΔP and generally is independent of other factors. In other words, a given ΔP corresponds to a certain removal efficiency and defines the relationship between cumulative contaminant loading, M, and ΔP, pressure drop, as

ΔP=A+Bexp(CM)

where:

-   -   A=constant and is assumed to be equal to the initial ΔP when         M=0.     -   B=constant that determines the point in time or mass loaded         where the ΔP transitions into exponential increase. The larger B         is, the sooner the curve goes exponential. B=f(concentration)         and the higher the concentration, the larger B is.     -   C=constant and reflects how steeply ΔP climbs during exponential         phase. The larger C is, the steeper the climb. With the y-axis         plotted on log scale, “C” indicates the slope of the line after         the transition point is passed, as determined by “B”.         C=f(concentration, type of contaminants) and the more important         soft contaminants are and the higher the concentration of soft         contaminants, the larger C is. When filter test data is used to         determine constants, this term can correct for differences in         contaminant concentration and type (hard vs. soft).

FIG. 7 illustrates the use of this relationship as an algorithm in a system for monitoring and indicating filter life as contemplated herein including automatic correction steps. Other alternative equation forms can be identified that closely match the observed ΔP versus time response (or the mass-loaded or fluid filtered response). For example, ΔP=A+B*M^(C) generates very similar response (with very different values of B, and C vs. the exponential relationship (see FIG. 5)); or a modified form of this equation, such as ΔP=A+B*exp(−M/C) which is functionally identical to the original form.

Alternatively, the pressure drop relationship can be expressed as a function of V, the volume of fluid filtered, as

ΔP=A+Bexp(DV)

D=Cc _(ave)

where: D=constant and reflects how steeply ΔP climbs during the exponential phase; and c_(ave)=average upstream contaminant concentration in the volume of fluid filtered. FIG. 8 illustrates the use of this relationship as an algorithm in a system for monitoring and indicating filter life as contemplated herein including automatic correction steps.

Because ΔP is a function of flow rate and fluid viscosity, it should be understood that, in order to use the previous equations the pressure drop data obtained for an operating system, pressure drop data may be normalized with respect to flow rate and viscosity. For convenience, it may be assumed in this description that the results may be normalized to the flow rate and viscosity conditions that were used under reference conditions, such as a laboratory, to generate the initial values for A, B, C and D. Moreover, it may be assumed that the ΔP varies proportionally in response to flow-rate and viscosity, as is typically found to be the case for viscous-dominated flow of liquids through filter devices where the filter element itself causes a majority of pressure drop. This assumption may be violated for poorly-designed filter housings containing “orifice-like” passages or constrictions that cause high inertial losses (dominating the total flow restriction), but those cases are not common. The normalized pressure drop, ΔP_(N), may be defined as

${\Delta \; P_{N}} = {\Delta \; P\frac{Q_{T}\eta_{T}}{Q\; \eta}}$

where:

Q_(T)=flow rate at reference condition

Q=flow rate

η_(T)=fluid viscosity at reference condition

η=fluid viscosity

In FIGS. 7 and 8, this relationship is used to normalize pressure drop data in order to apply the appropriate algorithms to calculate remaining useful life. In cases where the ΔP does not vary proportionally in response to flow-rate and viscosity, other appropriate normalization relationship may be used.

FIGS. 2-4 show how ΔP changes as a function of M, in response to adjustments in A, B and C. M can be related to remaining useful life as follows:

$M = {\frac{1}{C}{\ln \left( \frac{{\Delta \; P} - A}{B} \right)}}$ M_(R) = M_(T) − M $\frac{M_{R}}{M_{T}} = {\frac{t_{R}}{t_{T}} = {\frac{M_{T} - M}{M_{T}} = {1 - \frac{M}{M_{T}}}}}$ $t_{R} = {t_{T}\left( {1 - \frac{M}{M_{T}}} \right)}$

where:

-   -   M=cumulative contaminant loading     -   M_(R)=remaining contaminant holding capacity of filter     -   M_(T)=total contaminant loading capacity of filter     -   Q_(T)=flow rate used to determine capacity of filter     -   t_(T)=time required for filter to be fully loaded at Q     -   t_(R)=estimated time remaining before filter is fully loaded         with contaminant.

This can further be related to time, flow rate and volume filtered by:

M=Vc=Σ ₀ ¹ QcΔt

where:

-   -   V=Volume filtered     -   t=Time     -   Q=Flow rate. Q=f(t) and depends on duty cycle.     -   c=Concentration. c=f(t) and depends on the operating         environment, duty cycle, ingression, vibration, flow surges, and         other application factors.         In order to calculate t_(R), M must be known. M can be estimated         using the equations presented above. M is directly proportional         to V, assuming constant concentration. M can be estimated more         accurately if Q and/or c as a function of time can be determined         or is known. FIG. 1 provides a schematic diagram of         contamination sources that may provide bases for the algorithm.

Alternatively, t_(r) can be calculated from V as follows:

$V = {\frac{1}{D}{\ln \left( \frac{{\Delta \; P} - A}{B} \right)}}$ V_(R) = V_(T) − V $\frac{V_{R}}{V_{T}} = {\frac{t_{R}}{t_{T}} = {\frac{V_{T} - V}{V_{T}} = {1 - \frac{V}{V_{T}}}}}$ $t_{R} = {t_{T}\left( {1 - \frac{V}{V_{T}}} \right)}$

where:

-   -   V_(R)=remaining volume of fluid that can be filtered before the         filter reaches its contaminant holding capacity of filter     -   V_(T)=total volume of fluid that can be filtered by a new filter         before its contaminant loading capacity is reached     -   Q_(T)=flow rate used to determine capacity of filter time         required for filter to be fully loaded at Q     -   t_(R)=estimated time remaining before filter is fully loaded         with contaminant.

This can further be related to time, flow rate and volume filtered by:

where:

V=Σ _(u) ^(t) QΔt=Q _(cum) t

-   -   t=Time     -   Q=Flow rate. Q=f(t) and depends on duty cycle.     -   Q_(cum)=Sum of the flow rate values for each time increment         sampled by the system during current filter service interval.

Various corrections to the disclosed algorithm may be made over time in order to better estimate the remaining useful life of a filter make the algorithm more accurate with respect to estimating useful life. For example, a volume correction to the estimated remaining useful life may be performed based on deviation from expected average flow rate. A concentration correction to the estimated remaining useful life may be performed based on ΔP or knowledge of application. In a first instance, the algorithm assumes that removal efficiency is a function of ΔP as previously noted. An engine sensitivity correction to the estimated remaining useful life may be performed based on duty cycle. Duty cycle is a term used to describe the severity of the engine's operating conditions and application. It can be defined quantitatively or qualitatively in various ways. For example, quantitatively, duty cycle can be defined as the ratio of the time-weighted average produced power to the engine rated power; the percentage of time that the engine operates at rated power (or some fraction of rated power); the ratio of the time-weighted average fuel consumption rate to the fuel consumption rate at rated power output; the ratio of the time-weighted average exhaust temperature to the rated exhaust temperature at standard temperature and pressure (STP) ambient conditions; the ratio of the time-weighted turbocharger boost pressure or temperature to the boost pressure or temperature at rated power output conditions; the ratio of some function of time and temperature of oil sump or coolant to the same function operating at rated power conditions; or some combinations of two or more of the previously listed definitions. Duty cycle can also be expressed in descriptive terms, e.g., severe, moderate, light, using the previously mentioned quantitative definitions, or others, such as fuel economy. Duty cycle affects, the magnitude and/or frequency of flow transients, as well as temperature, which may be utilized to perform an empirical correction. Duty cycle affects the overall generation of contaminants, both soft and hard.

FIG. 6 shows an algorithm, in flow chart form, when the various constants of the algorithm are corrected manually, such as by service personnel. Examples of indications that t_(T) is too low and that constants B and/or C needs to be manually readjusted, include, but are not limited to: high ΔP at change-out; loss of power (fuel); smoke (fuel); poor oil analysis; collapsed pleats; a slimy or coated filter; unscheduled repairs; fuel economy low or other indicator of harsher duty cycle; operating in dirty/dusty/off-road environment; and fluid supply/suppliers known to have poor quality.

Examples of indications that t_(T) is too high and that constants B and/or C need to be manually readjusted, include, but are not limited to: ΔP low at change-out; Components in good shape at scheduled service interval; Very good oil analysis; Fuel economy high or other indicator of light duty cycle; Operating in clean environment; and Fluid supply/suppliers known to have high quality.

-   -   For the automatic version, illustrated in FIGS. 7 and 8 in flow         chart form, adjustments in constant values may be performed         automatically by the ECM or controller. Both service and         automatic algorithms can be further simplified, if desired, by         setting C or D as a true constant, A is a constant that is         approximately equal to the initial ΔP when M=0;     -   M is cumulative contaminant loading;     -   B is a constant that determines the point where the ΔP         transitions to an exponential rate of increase;     -   D=Cc_(ave)     -   C is a constant that reflects the exponential plugging slope;         and

c_(ave)=average upstream contaminant concentration in the volume of fluid filtered and adjusting only B. In this case, B is a correction for both concentration and type of contaminant. This would simplify the calculation with minimal loss of accuracy and make adjustment easier, especially for manual application.

For the automatic version using filtered volume calculation, illustrated in FIG. 8 in flow chart form, adjustments in constant values may be performed automatically by the ECM or controller and the filter serviced as summarized. As shown in block 1, the first time when the system is used, initial values for A, B, D, ΔP_(T) (the terminal drop pressure), Q_(T), η_(T) based on reference conditions are entered into the system. The terminal pressure drop, ΔP_(T), is the pressure drop at which the filter should be replaced or serviced. In FIG. 8, blocks 2-5 indicate steps performed when a new filter is installed in which initial conditions for the new filter are entered and constant A is set based on the measured pressure drop, flow rate and viscosity. Blocks 6-10 and diamond 8 show measurement steps and calculation of the remaining volumetric capacity of the filter, V_(R). Blocks 11-14 and diamonds 12-13 show the portions of the algorithm that determine whether or not the filter should be replaced or serviced. Diamonds 15-16 show the portions of the algorithm that determine whether or not the constants B and/or D need to be recalculated. Block 17 shows the portion of the algorithm in which the constant B is recalculated. Block 19 shows the portion of the algorithm in which the constant D is recalculated. Diamonds 18 and 20 show the portions of the algorithm that determine whether or not the filter is damaged and needs to be replaced. Other means of implementing the algorithms to optimize filter service intervals are anticipated, such as by storing historical values for A, B, C and/or D and comparing them to present values in order to determine whether or not the filter has been damaged or compromised or otherwise be replaced.

In some embodiments, the disclosed systems and methods may include or utilize hardware components as illustrated in FIG. 9. Specifically, the systems and methods may include or utilize one or more components selected from: a filter; a ΔP sensor to measure ΔP across the filter; a flow sensor to measure or estimate the flow rate through the filter; a timer to measure time at given flow rate or engine condition; a viscosity and/or temperature sensor to measure or estimate the viscosity of the fluid; a controller to accept sensor output, perform the calculations required by the algorithm, and output a signal indicative of the results of the algorithm calculation; an appropriate output device or means; and a sensor or other means to sense pressure flow rate, viscosity and/or temperature data which may be inputted into a device that utilizes the algorithm disclosed herein to correct pressure drop. Suitable sensors may include, but are not limited to, a flow meter, a separate ΔP sensor; a temperature sensor; a viscosity sensor; and a virtual sensor that makes the correction based on engine operating conditions.

In the foregoing description, certain terms have been used for brevity, clearness, and understanding. No unnecessary limitations are to be inferred therefrom beyond the requirement of the prior art because such terms are used for descriptive purposes and are intended to be broadly construed. The different configurations, systems, and method steps described herein may be used alone or in combination with other configurations, systems and method steps. It is to be expected that various equivalents, alternatives and modifications are possible within the scope of the appended claims. It is provisionally and preliminarily expected that the claims will include at least the following. 

1. A method for determining and indicating remaining useful life of a filter for filtering fluid in a given machine application, comprising utilizing an algorithm for calculating remaining useful filter life, providing initial constants for the algorithm based on at least one of laboratory performance and historical performance, adjusting the constants during filter usage based on at least one of actual performance of the filter and behavior of the machine.
 2. The method according to claim 1, wherein the method optimizes cost-efficiency and value to a customer by maximizing filter service intervals and minimizing available but unused filtration capacity prior to end-of-filter-life.
 3. The method according to claim 1, wherein the method maintains fuel economy by reducing or limiting total work required to pump fluid through the filter.
 4. The method according to claim 1, wherein the algorithm is utilized to calculate remaining useful filter life based on ΔP.
 5. The method according to claim 1, wherein the algorithm is utilized to calculate remaining useful filter life based on filter efficiency.
 6. The method according to claim 1, wherein the algorithm is utilized to calculate remaining useful filter life based on particle concentration downstream of the filter.
 7. The method according to claim 1 wherein the machine is an internal combustion engine.
 8. The method according to claim 1, wherein the machine is a hydraulic machine.
 9. The method according to claim 1 wherein the engine has an ECM, electronic control module, and wherein the constants are adjusted during operation based on data obtained or processed by the ECM including historical behavior data.
 10. The method according to claim 9 comprising using ECM data to calculate the volume of fluid flow through the filter, using the algorithm to indicate remaining useful filter service life, and determining whether filter service is required.
 11. The method according to claim 1, wherein the constants are adjusted during operation based on data obtained by a control module that is separate from the machine including historical behavior data, which data optionally is transmitted to a remote computer or controller that utilizes the algorithm.
 12. The method according to claim 1 comprising manually adjusting the filter constants.
 13. The method according to claim 12, wherein the filter constants are adjusted manually during operation.
 14. The method according to claim 12, wherein the filter constants are adjusted manually at a service center.
 15. The method according to claim 1 comprising automatically adjusting the constants during operation.
 16. The method according to claim 1 comprising adaptively varying and adjusting the constants during operation.
 17. The method according to claim 16, comprising adjusting the algorithm according to at least one of contaminant concentration, contaminant type, duty cycle, and filter type.
 18. The method according to claim 1 comprising observing performance of the filter during operation, and adjusting the constants during operation based on observed filter performance.
 19. The method according to claim 1 comprising observing performance of the filter during operation, and adaptively updating the algorithm either manually or automatically based on observed filter performance.
 20. The method according to claim 1 comprising providing an indication alert to a user that the filter is approaching an end-of-useful-life condition based on total volume of fluid filter and based on operating conditions to which the filter has been exposed.
 21. The method according to claim 1 wherein the filter is a fuel water separator.
 22. The method according to claim 1 comprising adjusting the constants during operation according to a number of start-up and shut-down events of the machine.
 23. The method according to claim 1 comprising using the algorithm as an adaptive tool for scheduling filter changes so that a filter is not changed too soon, thereby increasing cost of operation, nor too late, thereby resulting possible damage to the machine.
 24. The method according to claim 1 comprising using the algorithm as an adaptive tool which can be adjusted by service personnel or automatically adjusted based on field experience to more accurately estimate remaining useful life of the filter.
 25. The method according to claim 1 comprising determining the volume of fluid filtered by the filter, and using the algorithm to calculate remaining useful life of the filter based on fluid volume filtered.
 26. The method according to claim 25 comprising performing the determining step by calculating the volume of fluid filtered by the filter by tracking flow rate and time.
 27. The method of claim 1, wherein the algorithm has a formula selected from: ΔP=A+Bexp(CM), ΔP=A+B*M ^(C), and ΔP=A+B*exp(−M/C) wherein: ΔP is pressure drop; M is cumulative contaminant loading; A is a constant that is approximately equal to the initial ΔP when M=0; B is a constant that determines the point where the ΔP transitions to an exponential rate of increase; and C is a constant that reflects the exponential plugging slope.
 28. The method of claim 1, wherein the algorithm has a formula: ΔP=A+Bexp(DV) wherein: A is a constant that is approximately equal to the initial ΔP when M=0; M is cumulative contaminant loading; B is a constant that determines the point where the ΔP transitions to an exponential rate of increase; D=Cc_(ave); C is a constant that reflects the exponential plugging slope; and c_(ave)=average upstream contaminant concentration in the volume of fluid filtered.
 29. A system for determining and indicating useful life of a filter for filtering fluid in a given machine the system using an algorithm for calculating remaining useful filter life and comprising: (a) a filter; (b) a ΔP sensor to measure ΔP across the filter; (c) a controller to accept sensor output, perform the calculations required by the algorithm, and output a signal indicative of the results of the algorithm calculation; (d) an output device; and (e) one or more sensors to correct pressure drop data for changes in viscosity and flow; wherein, in the system, initial constants for the algorithm are provided based on at least one of laboratory performance and historical performance, and the initial constants are adjusted during filter usage based on at least one of actual performance of the filter and behavior of the machine.
 30. The system of claim 29, wherein the algorithm has a formula selected from: ΔP=A+Bexp(CM), ΔP=A+B*M ^(C), and ΔP=A+B*exp(−M/C) wherein: ΔP is pressure drop; M is cumulative contaminant loading; A is a constant that is approximately equal to the initial ΔP when M=0; B is a constant that determines the point where the ΔP transitions to an exponential rate of increase; and C is a constant that reflects the exponential plugging slope.
 31. The system of claim 29, wherein the algorithm has a formula: ΔP=A+Bexp(DV) wherein: A is a constant that is approximately equal to the initial ΔP when M=0; M is cumulative contaminant loading; B is a constant that determines the point where the ΔP transitions to an exponential rate of increase; D=Cc_(ave); C is a constant that reflects the exponential plugging slope; and c_(ave)=average upstream contaminant concentration in the volume of fluid filtered. 